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Abstract 

The high-lights of ground-based very-high-energy (VHE, E > 100 GeV) gamma- 
ray astronomy are reviewed. The summary covers both Galactic and extra- 
galactic sources. Implications for our understanding of the non-thermal Uni- 
verse are discussed. Identified VHE sources include various types of supernova 
remnants (shell-type, mixed morphology, composite) including pulsar wind neb- 
ulae, and X-ray binary systems. A diverse population of VHE-emitting Galactic 
sources include regions of active star formation (young stellar associations ), 
and massive molecular clouds. Different types of active galactic nuclei have 
been found to emit VHE gamma-rays: besides predominantly Blazar-type ob- 
jects, a radio-galaxy and a flat-spectrum radio-quasar have been discovered. 
Finally, many (presumably Galactic) sources have no convincing counterpart 
and remain at this point unidentified. A total of at least 70 sources are currently 
known. The next generation of ground based gamma-ray instruments aims to 
cover the entire accessible energy range from as low as ~ 10 GeV up to 10 5 
GeV and to improve the sensitivity by an order of magnitude in comparison with 
current instruments. 



1 Introduction 

The highest energy photons known are produced in astrophysical processes involv- 
ing even more energetic particles presumably accelerated through stochastic accel- 
eration mechanisms as suggested initially by |Fermi| ([T949 ). Therefore, observations 
of VHE photons provide a direct view of the astrophysical accelerators of charged 
particles and allow to identify the individual sources of cosmic rays: VHE photons 
open our view to the "accelerator sky". The origin of the Galactic population of 
charged cosmic rays has remained since its discovery by V. Hess in 1912 until today 
a long-standing question of astro- and particle physics. A widely favored model on 
the origin of cosmic rays assumes that diffusive shock acceleration takes place in the 
expanding blast waves of supernova remnants converting 10 — 20 % of the kinetic en- 
ergy into cosmic rays. Under this assumption, Galactic supernova remnants provide 
sufficient power (i=s 10 41 ergs s _1 ) in order to balance the escape losses of Galactic 
cosmic rays as well as produce a power-law type distribution of particle energy that 
closely resembles the cosmic ray spectrum measured locally (see e.g. Ginzburg & 
Syrovatskii| 1964] |Hillas|2006 1. 



1 



Observations of VHE-gamma-rays, mainly by imaging air Cherenkov telescopes in 
the last decade, have surpassed the anticipated detection of a few supernova rem- 
nants and have established a rich and diverse collection of VHE sources. Specifi- 
cally, the currently active generation of imaging air Cherenkov telescopes (H.E.S.S. 1 , 
MAGIC 2 , and VERITAS 3 ) have fulfilled and by far exceeded the expectations that 
were based upon the pioneering previous generation of experiments: the results ob- 
tained in the last years have shown that VHE emission is common to a variety of 
different source types - not only shell-type SNRs. 

The current view of the source distribution in the VHE sky is shown in Figure [T] 
where all known sources are displayed (status of September 2008). Note, the sen- 
sitivity achieved varies greatly across the sky. The best sensitivity is reached in the 
inner Galactic disk where a dedicated survey with the air Cherenkov telescopes of 
the H.E.S.S. experiment has been performed (see section 2.2 1. 
The most remarkable feature of the gamma-ray sky is not evident in this picture: 
each source shown is an accelerator of particles up to and beyond TeV energies ("ac- 
celerator sky"). A similar conclusion can not be drawn when looking at source pop- 
ulations detected in other wavelength bands. This makes the VHE band a sensitive 
window to detect non-thermal particle accelerators. Future neutrino telescopes will 
almost certainly have the potential to detect some of the brighter VHE sources (and 
discover new sources that are optically thick for VHE gamma-rays). The detection 
of high energy neutrinos is experimentally a challenging task and requires tremen- 
dous efforts (see e.g. the Ice-CUBE neutrino telescope in the Antarctic ice). How- 
ever, the observation of a neutrino source will ultimately demonstrate the presence of 
accelerated nuclei - largely model-independent. Different to the clear observational 
signature in the neutrino channel, VHE gamma-rays are sensitive to accelerated elec- 
trons (through inverse Compton scattering) as well as to accelerated nuclei (through 
neutral meson production and decay). 

The scope of this article is limited to observational highlights obtained with 
ground based instruments and does not aim at presenting a complete review of the 
field of VHE astrophysics. The paper is structured in the following way: In sec- 
tion |2j the experimental techniques of ground-based instruments (both imaging and 
non-imaging) are described before continuing with the census of today's VHE sky in 
section|3] Section|4]provides a a short overview on VHE gamma-rays as probes of the 
interstellar medium including Lorenz-invariance violating effects related to structure 
of space-time at Planck-scales. Finally, in section [5] this review is concluded with 
some comments on the future of the field. 



2 Observational techniques 

Non-thermal gamma-ray sources produce typically power-law type energy spectra in 
which the flux drops with increasing energy. As a consequence of this, space based 
detection techniques are limited by the small detection rate of photons at high ener- 



1 http://www.mpi-hd.mpg.de/hfin/HESS/HESS.html 

2 http ://w w wmagic . mppmu.mpg.de/ 



3 http://veritas. sao.arizona.edu/ 
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Figure 1 : All VHE sources as known in September 2008 are shown as circles in a 
Hammer- Aitoff projection of a Galactic coordinate system. The Galactic center is in 
the center of the picture. 
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gies (E < 10-100 GeV). The detection of VHE (E > 100 GeV) photons requires 
collection areas that are substantially larger than the ff(m 2 ) of satellite based de- 
tectors. One approach towards large collection areas are ground based observation 
techniques which however are limited to energies sufficiently high (E > 5 GeV) to 
produce a detectable air shower. 

Ground based gamma-ray observations are always based upon air shower detection 
techniques. High energy particles (photons, electrons, and nuclei) initiate exten- 
sive air showers in the atmosphere which acts as a natural calorimeter: secondary 
particles are detected either when they reach the ground or by the Cherenkov light 
produced in the atmosphere. At ultra-high energies, it becomes feasible to detect 
fluorescence light as well as radio signals which are not discussed here (see e.g. the 
contribution by J. Horandel at this conference). 

Generally, the air shower detection technique relies on one of the following ap- 
proaches: 

• Shower-front sampling of particles 

• Shower-front sampling of Cherenkov light 

• Imaging of the air shower using Cherenkov light. 

The benefits (marked '+') and draw-backs ('-') of the shower-front sampling of par- 
ticles and the imaging air Cherenkov technique are mainly the following 4 : 

• Imaging technique: 

+ sensitivity (rs 10~ 13 ergs/(cm 2 s) _1 at TeV-energies, 50 hrs exposure) 
+ angular resolution (few arc minutes per event) 
+ spectroscopy (AE/E w 15 %) 
+ low energy threshold 

- field of view (a few degrees, w msrad) 

- duty cycle (w 10 %) 

• Shower front sampling of particles: 

+ high duty cycle (w 95 %) 
+ large field of view (w 2 srad) 

- high energy threshold 

- sensitivity (w 10~ 12 ergs cm~ 2 s _1 , 1 year exposure) 

- angular resolution (0.3 — 0.7°) 

- energy resolution (AE/E w 30 — 100 %) 

4 Shower- front sampling of Cherenkov light is discussed seperately below. 
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The obvious complementarity of the two techniques justifies the further development 
of both of them in the future: non-imaging techniques cover simultaneously a large 
field of view and can be operated 24 hours a day. On the other hand, air Cherenkov 
telescopes have a narrow field of view (a few milli steradians), reaching a supreme 
sensitivity, but can only be operated in clear, dark nights. 

In addition to the indirect ground based techniques discussed until now, high en- 
ergy gamma-rays can be detected above the atmosphere directly by pair-conversion 
detectors. With these detectors, events initiated by charged particles can be effec- 
tively suppressed by e.g. a scintillation veto shield. Obviously, this is not possible 
for ground based indirect detection techniques where a separation of gamma-ray 
and cosmic-ray events has to be done using the information of the air shower itself 

(gamma-hadron separation). 

With the commissioning of the Fermi satellite 5 ( Atwood et al.|2007 1, ground-based 
gamma-ray astronomy will benefit from the simultaneous operation of essentially a 
sensitive all-sky monitor operating at lower energies. The sensitivity of Fermi is well 
matched by the sensitivity of ground-based experiments so that the two instruments 
will be able to provide detection or constraints on the spectral shape across 6-7 orders 
of magnitude in energy. 



2.1 Shower-front sampling techniques 

Particle arrays. The experimental approach for the measurement of cosmic rays as 
well as the search for gamma-ray sources had been largely dominated by air shower 
arrays until the beginning of the 1990's. These arrays of particle detectors measure 
the arrival time of the shower front and the lateral particle density distribution. The 
total number of particles in the air shower is used to reconstruct the total energy of 
the primary particle. The relative timing of the individual particle detectors in the 
array is used to determine the direction of the air shower (typical angular resolutions 
of 1° can be reached) while both the timing and particle density measured on the 
ground are in principle useful to discern electromagnetic from hadronic air-showers 
(gamma-hadron separation). However, the gamma-hadron separation is naturally 
limited by the filling factor 6 of the detector array which typically amounts to less 
than 1%. An important development in the field was the innovation introduced by 
the MILAGRO collaboration (see e.g. Atk ins et al.||200 4): they use a pond of wa- 
ter which is instrumented with photo-multiplier tubes (PMTs). The PMTs detect 
Cherenkov light from the secondary particles entering the pond's water. This way, 
the filling factor is increased to almost 100% and the entire particle distribution is 
sampled by the detector. Using the dumpiness of the particle distribution as an indi- 
cator for a hadronic air shower 7 , it has been possible to reach sufficient sensitivities 
to detect VHE gamma-ray sources. 

Future projects along this direction aim at installing such a detector at high altitude, 
where the energy threshold can be decreased to reach values below 1 TeV. Parallel 

5 formerly known as GLAST 

6 the fraction of the total surface covered with active detector surface 

7 for electromagnetic air showers, the lateral distribution is smoother than for a hadronic shower which 
contains muons and sub-showers 



5 



to the MILAGRO group, the ARGO collaboration has developed a new type of air 
shower detector using resistive plate chambers to increase the filling factor. Their 



installation is already located at high altitude and is starting operation ( Aielli & The 



|Argo-YBJ Collaboration 2007 1. For a review of ground based non-imaging detec- 
tors, see [Lore nz (2007}. 



Shower-front sampling with Cherenkov light. The experimental technique used 
by non-imaging air Cherenkov detectors marks the transition from non-imaging to 
imaging observations. These detectors sample the arrival time and density of air 
Cherenkov photons using either open PMTs like e.g. the THEMISTOCLE exper- 
iment ( |Fontaine et~aLl|1990| l, AIROBICC array ( |Karle et alT]|1995) > or large reflect- 
ing surfaces as e.g. heliostats provided by solar power plants. The latter approach 
has been pursued by a number of groups including the C.A.C.TU.S. (|Marleau et aL 



|2005) , CELESTE ( |Pare etal.|2002| as well as STACEE ( |D. M. Gingrich et al.|2 005 ). 
While arrays of open PMTs retain the large field of view of classical air shower sam- 
pling arrays, the heliostat arrays have a very small field of view but a substantially 
smaller energy threshold reaching below 100 GeV as achieved e.g. by the CELESTE 
experiment . The potential for discovering faint sources is limited by the gamma- 
hadron rejection power of these shower sampling instruments. 



2.2 Imaging air Cherenkov telescopes 

The field of ground-based gamma astronomy has been largely driven by the remark- 
able results obtained with the imaging air Cherenkov telescopes (IACTs). Extensive 
air showers emit in the forward-direction a beam of atmospheric Cherenkov light 
with an opening angle of « 1°. This beam illuminates almost homogeneously a cir- 
cular region on the ground with a diameter of 200-300 m (depending on the altitude 
and inclination of the shower axis). An optical telescope pointing parallel to the 
shower axis and located within the illuminated footprint of the shower can make an 
image of the air shower against the background light of the night sky, provided the 
camera is sufficiently fast to integrate the short Cherenkov flash of only a few nano- 
seconds. The image provides information on the original particle's energy, direction, 
and on its nature (nucleus or photon). The reconstruction of these parameters is im- 
proved considerably when more than one telescope is used in a stereoscopic set-up 
where the telescopes are separated by 50-150 m in order to provide a baseline for 
triangulating the atmospheric air shower. The stereoscopic technique has become 
the nominal standard for all current and future installations. 

The performance of today's telescopes is essentially characterized by the sensitiv- 
ity to detect VHE sources with an energy flux down to 10~ 13 ergs (cm 2 s)~' in 
50 hrs of observation time; this corresponds to a minimum detectable luminosity of 
^min ~ 10 31 ergs s _1 (d/1 kpc) 2 at a distance of 1 kpc or more suitable for extra- 
galactic objects L min w 10 41 ergs s -1 (c//100 Mpc) 2 at a distance of 100 Mpc. 
The angular resolution of each reconstructed primary y-ray is typically better than 
6 arc min. The relative energy resolution is comparably good and reaches values of 
AE/E w 10 — 20 %. This is sufficient to detect and characterize spectral features like 
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curvature or even lines from e.g. self-annihilation of Dark matter particles. For a list 
of currently operating IACTs see e.g. |Hinton| ( |2007| ). 

3 A census of today's VHE gamma-ray sky 

The currently known list of VHE sources encompasses more than 70 sources (see 
Fig . [TJ . The number of sources is modest in comparison with catalogues assembled 
at lower energies. However, it is remarkable, how many different types of objects 
are actually emitting VHE gamma-rays and are therefore accelerators of multi-TeV 
particles. 

Before discussing the different source types that have been detected, it is enlighten- 
ing to list promising source types which have not been detected in the VHE band so 
far (the possible reason(s) for non-detection listed in parenthesis): 

• Star burst galaxies (too faint): The high star forming rate should lead ulti- 
mately to an enhanced rate of supernova explosions. The upper limits on VHE 
emission from e.g. NGC 253 ( Aharonian et al. 2005h) are still higher than the 
expected flux (Domingo-Santamarfa & Torres|20 05 ) but deeper observations 
(ca. 50 hrs) will eventually reach the required sensitivity to either detect a 
signal or put meaningful constraints on the models 8 . 

• Gamma-ray bursts (too short, too far): These extremely powerful explo- 
sions are transient events requiring fast turn-around times and/or wide-field- 
of-view instruments in order to capture these objects during their outbursts. 
No convincing GRB detection has been claimed so far. 9 The ultra-fast slew- 
ing MAGIC telescope has succeeded in observing (but not detecting) a GRB 
while the outburst was still ongoing ( |Albert et al.|2006c| ). Non-imaging all-sky 
survey instruments like MILAGRO are probably more likely to detect these 
transient events including also the distinct class of short GRBs (for constraints 



see e.g. Abdo et al. 2007c i. The visibility of distant GRB events is limited due 
to absorption of the energetic photons emitted at cosmological distances in 
pair-production processes with the optical-to-infra-red background light. An 
energy threshold well below 100 GeV is crucial to be able to observe GRBs at 
red-shifts z > 0. 1 . 

Clusters of galaxies (too extended, too faint): So far, no group or cluster 
of galaxies has been detected in the VHE band. Clusters of galaxies confine 
effectively cosmic rays and are expected to produce VHE-emission via inelas- 
tic scattering of nuclei with the intra-cluster gas (Volk et al. 1996, Pfrommer 



& EnBlin 2004). In addition to cosmic rays injected by normal galaxies and 
accelerated in large-scale shocks, AGN activity could also lead to an addi- 



tional contribution to the intra-cluster cosmic ray density ( Hinton et al. 2007b ). 



8 The CANGAROO collaboration published a claim for a detection from the starburst galaxy NGC 253 
which was however later retracted jltoh et al.|2007| . 

'Some evidence for GRB detections have been published in the past (see e.g. |Padilla et al.|l998[|Atkins| 
|et al.|2000} . 
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Current upper limits do not constrain severely the non-thermal population of 
cosmic rays in the clusters (Dom ainko et al.|2007) >. 

• Pulsars (early cut-off in energy spectrum): Fast rotating neutron stars are 
expected to give rise to acceleration within the magnetosphere (for a review, 
see e.g. |Kaspi et"aL||2006| ). This leads to pulsed emission (mainly from cur- 
vature radiation) that has been detected at least from six isolated pulsars up to 
GeV energies with the EGRET spark chamber detector on board the Comp- 
ton Gamma-Ray Observatory (Thompson et al. 1999). However, the energy 
spectra of pulsars are also expected to show a sudden cut-off at a few GeV 
which has made them so far invisible to ground-based telescopes with an 
energy threshold around 100 GeV. Most searches for pulsed emission above 
100 GeV have so-far not been successful in finding VHE emission from pul- 
sars - consistent with the expectations ( |Aharonian et al.||2007g| | Albert et al.| 
2007b). The MAGIC collaboration has recently reported the first detection of 
pulsed emission with ground based instruments above 25 GeV at the level of 
6.4 a (Teshima 2008). Additionally, pulsed emission has been suggested to 



be produced via bulk Compton-scattering processes in the un-shocked pulsar 
wind-zone (Bogovalov & Aharonian 2000). Non-detection of pulsed emission 
from the Crab pulsar has been used to constrain the formation region of the 
ultra-relativistic wind (Aharon ian et al.|20 04a). 

The breakdown of known VHE sources includes 50 Galactic objects: 5 shell-type 
SNRs, 2 mixed-morphology SNRs, 2 composite SNRs, 20 pulsar-wind nebulae, 2 
stellar associations, 4 X-ray binary systems, and roughly 18 sources without clear 
association to known objects. The remaining 23 objects are of extra-galactic origin: 
21 Blazars, one Fanaroff-Riley Type I (M 87), and one flat-spectrum radio quasar 
(3C279). 



3.1 Shell-type supernova remnants 

Shell-type supernova remnants are commonly considered to be the best candidates 
to accelerate Galactic cosmic rays. During the supernova explosion (either a thermo- 
nuclear explosion or a core-collapse event), the stellar atmosphere is ejected with an 
initial velocity of a few 1000 km s _1 carrying 10 51 erg in kinetic energy and driv- 
ing a shock front in the interstellar medium. While slowing down during the Sedov 
phase 10 , the expanding shock front heats up the ambient medium giving rise to ther- 
mal X-ray emission. The shock is expected to dissipate a good fraction (10 — 30 %) 
of its energy during the Sedov phase in the form of accelerated particles (electrons 
and nuclei). In this case, the total power injected by roughly 1 supernova explosion 
every 30 years is sufficient to balance the escape losses of Galactic cosmic rays. 
Radio- as well as X-ray synchrotron emission is detected and attributed to electrons 
accelerated by the forward shock which in projection produces the characteristic 
shell-like morphology. For a few objects, the non-thermal X-ray component domi- 
nates entirely the observed X-ray spectrum (e.g. SN 1006, RX J1713.7-3946). 

I0 once the SNR has swept up a comparable amount of matter from the interstellar medium as its ejecta 
mass, the free expansion phase ends and the Sedov phase of SNR evolution starts 
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Figure 2: Upper panel, from left to right: RX J0852.0-4622 ("Vela Jr") ( Aharonian| 
etaL|2007f), RX J1713.7-3946 (|Aharonian et~aL|2006d), bottom panel: RCW 86 



( Hoppe et a 



^2007) , SN1006 (Naumann-Godo & H.E.S.S. CoUaboration|2008) . The 



color-coded images are the excess maps as obtained with the H.E.S.S. telescopes (see 
text for references). Overlaid in white contours (RX J0852. 0-4622) is the smoothed 
ROSAT image obtained above 1.3 keV. The black contours overlaid on the VHE- 
image of RX J1713.7-3946 follow the ASCA image (1-3 keV). The white contours 
overlaid on RCW 86 are the significance contours (starting at 3 standard deviations, 
increasing by 1) The white contours on the preliminary excess map from SN 1006 
follows the Chandra X-ray map smoothed to match the H.E.S.S. point-spread func- 
tion. 
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However, radio and X-ray observations are sensitive mainly to electrons and only in- 
directly to the presence of energetic nuclei: X-ray observations with sub-arcsecond 
spatial resolution have revealed indirectly, that the forward shock in many SNRs is 
very likely modified by cosmic ray streaming instabilities leading to enhanced mag- 
netic field compression ( |Lucek & Bell|2 000 ). This in turn leads to observable effects: 
electrons produce X-ray synchrotron emission in small regions with increased mag- 
netic field upstream of the shock that appear as narrow X-ray filaments (d < 1 pc) 
that have been first detected using the superior angular resolution of the Chandra 
X-ray telescope (Bamba et al. 2003). The inferred magnetic field strength in the 



X-ray filaments of e.g. SN1006 is of the order of 100 £lG ( |Berezhko efaij|20"03" 



which is much larger than in the case of a linear shock thus indirectly revealing effi- 
cient acceleration of nuclei. Recent observations of fast variability of X-ray emission 
in filamentary structures seems to provide additional evidence for strong magnetic 
fields in RX J1713.7-3946 ( |Uchiyama et al.|2007| l. 

VHE gamma-ray observations are sensitive to both, electrons as well as nuclei. Suffi- 
ciently energetic electrons can radiate VHE gamma-rays through inverse Compton- 
up-scattering of soft seed photons (leptonic origin of VHE-emission) while nuclei 
produce neutral mesons decaying into gamma-rays in inelastic scattering events with 
the ambient gas (hadronic origin of VHE-emission). 

In the hadronic scenario, the energy loss-time for accelerated protons (T pp ^ K o ~ 
4.5 x 10 15 n- 1 s) varies only slowly with energy and depends mainly on the am- 
bient medium density (n# = n cm~ 3 ) which includes cold molecular and hot ion- 
ized gas as well as the ejecta of the progenitor star. The actual density of gas 
near the shock can be inferred e.g. by modeling thermal X-ray emission from 
heated and partially ionized gas. However, for SNR without thermal X-rays, the 
density can only be constrained rather loosely and uncertainties close to an or- 
der of magnitude can remain. The typical integrated luminosity in VHE gamma- 
rays (1-10 TeV) is from a population of energetic protons with total energy W p 
L y « W p X- 1 = 3.7 • 10 33 • 77/O.I • (£sn/10 51 ergs)« ergs s" 1 for an efficiency of 
shock acceleration 77 =0.1, ie. the fraction of the kinetic energy of the blast wave 
converted into charged particles accelerated 1 1 . A similar result on the SNR lumi- 
nosity has been obtained by |Drury et al.| ( [T994| l. 

Currently four VHE-emitting shell-type SNRs are detected and upper limits have 
been derived for a few more objects (see Table[T]for a list). Except for Cassiopeia A, 
the spatial extension of the SNR has been resolved (see Fig. |2j. Generally, the ob- 
served luminosities of the shell-type SNRs are consistent with the expected luminos- 
ity in the hadronic scenario. For RX J0852. 0-4622 (Vela Jr.), the distance estimates 
are controversial and therefore the luminosity uncertainties can be as large as 2 or- 
ders of magnitude. See the references given in Table[T]for more details. Within this 
uncertainty on the distance of Vela Jr., the VHE observations are consistent with the 
indirect evidence from X-ray observations for efficient shock acceleration of nuclei. 
A leptonic origin can not be excluded but is disfavored as it would lead to a com- 
parably small efficiency for acceleration of nuclei which in turn would not produce 
non-linear modifications as are observed in X-rays. For Cassiopeia A, the average 



1 The accelerated spectrum is assumed to follow a power-law with dn/dE oe E 
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magnetic field is known to be larger than 80 jUG|Cowsik & Sarkar (1980} and may 



be as large as 1 jj,G if magnetic field and particle energy density are in equipartition. 
This implies that the corresponding particle number density and the expected inverse 
Compton flux is too small to explain the observed VHE emission which effectively 
rules out a leptonic origin of the observed VHE emission in this particular case. 
A second crucial test of the SNR origin of cosmic rays is the spectral shape and 
maximum energy of particles accelerated 12 . The observed VHE energy spectra from 
RXJ17 13.7-3946 ( |Aharonian et al.|2006d| l and RX J0852.0-4622 jAharonian efaL] 



2007f ) indicate that the accelerated proton spectrum follows an E~ 2 type power-law 
and cuts off at energies of roughly 100-200 TeV which is about one order of magni- 
tude smaller than the energy of the "knee" in the cosmic ray all-particle spectrum at 
a few 10 15 eV. This discrepancy can be accommodated if these two SNRs are already 
too old and the shock has decelerated such that the maximum energy has dropped to 



the currently observed value (Ptuskin & Zirakashvili 2005) while the more energetic 



particles have already escaped the remnant. 
In this context, the young SNR Cassiopeia A 13 and its VHE properties are crucial to 
our understanding of the SNR contribution to cosmic rays. Cassiopeia A has been 
confirmed by the MAGIC and VERITAS collaborations to be a VHE-emitting source 
([Albert et al.|2007f||Humensky & VERITAAS Collaboration|2008|> after the initial 



discovery with HEGRA ( Aharonian et al.f2001a) . The available spectra measured 



with HEGRA and MAGIC are in agreement with each other but require a rather 
soft energy spectrum with a photon index close to 2.5. The total energy in protons 
W p »2x 10 49 ergs is only a small fraction of the kinetic energy of the blast wave and 



only a factor of 4-8 larger than the inferred energy in accelerated electrons ( Atoyan 



et al.||2000] >. Given the young age of the source which is at the beginning of the 



Sedov phase, the small value of W p may still be re-reconciled with Cassiopeia A to 
be a typical SNR. It will be quite interesting to measure with better accuracy the 
energy spectrum to higher energies in order to probe the maximum energy of ac- 
celerated protons which should be close to PeV energies. Furthermore, as indicated 
in Table [T[ future detections (or stronger constraints) of VHE-emission from other 
historical SNRs like Tycho or Kepler will provide important clues on cosmic -ray 
acceleration in these objects. 

3.2 Mixed morphology supernova remnants 

Besides young (f w ^(1000) yrs) shell-type supernova-remnants, VHE-gamma-ray 
emission has also been detected from the direction of mixed-morphology (MM) type 
supernova-remnants (SNRs) which are considered evolved systems. MM systems 
are characterized by a shell-like emission observed only in the radio-band with ther- 
mal, centrally-peaked X-ray emission predominantly from the interior of the SNR 
(see e.g Rho & Petre|1998) . The mechanism responsible for heating the gas is not 
well-understood. It has been suggested that mixed morphology SNR have entered a 
post-Sedov stage of their evolution which is characterized by an equal electron and 
ion temperature ( Kawasaki et al.|[2005 1. The fact that roughly 20 % of the known 



I2 which can not be inferred from observations in other wavelengths 
13 and also Vela Jr if the age is in the range of 200-300 yrs 
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Table 1: List of supernova-remnants (shell-type and mixed-morphology) observed 
in VHE gamma-rays . 



Name 


distance 




age 


Type 




[kpc] 


[10 3J ergs/s] 


[kyrs] 




Cas A 


3.4 fl 


l.4 b 


0.33 c 


Shell, II 


RX J 171 3.7-3946 




5J E 


1.6 e 


Shell, Wlbf 


VelaJr 


0.2* 


0.3° 


0.7* 


Shell, II 




0.33 /! 


0.6 G 


0.66^ 






1-2'' 


6-26 G 


4-19'' 




RCW 86 


2.8' 


5.5* 


1.8 (SN1857V 


Shell, ? 


SN 1006 


2.2' 


1.4 L 


1 


Shell, la 


W28 


1.9' 


0.5"' 


33 m 


MM 


IC443 


1.5" 


0.4° 


20^ 


MM 


Without detection (historical SNe) 


Tycho 


2.2* 


< 0.5'' 


0.4 


Shell, la 




4.5 s 


<2 r 


(SN1572) 




Kepler 


4.8" 


< 2.4 V 


0.4 


Shell, Ia w 




6.4" 


<4.2 V 


(SN1604) 





faeed et afl|l995l lAharonian et al.| 

,iirm q t*t 141 QQ7 l E Ahr,,™:™ ft al IJOi 



N2001aMAlbert et al.|2007f| , c Thorstensen et al. (20011. 
Koyama et al. ( 1997 1, £ Aharonian et al. l2006dl, e Koyama et al. ( 1997 1, f Cassam-Chenai et al. 12004), 
HAharonian et al.|j2007ft. ^lAschenbach et al.|jl999t, "iBamba et al.N2005t, 'I Slane et al. (2001 1 , J |Vink 
|et al.»2006), "|Hoppe et al |j2007), L |Naumann-God6 & H.E.S.S. Collaboration|j2008) ' |Velazquez et al 
<2002), "lAharonian et al.|j2008e), "|Claussen et al 1(1997), °|Albert et al.H2007c), P|Lee et al H2008) , 
^|Albinson et al. | |1986) , " |Aharonian et al »2001b^ "|Schwarz et al Hl995) , ^Winkler et al.||2003) , " 
|Reynoso & Goss 1 1999) , " | Aharonian eTaTf 2008b) '"' [Reynolds et al.[ 2007) 
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Figure 3: Left panel (a): The integrated CO emission in the velocity range 10- 
20 km/s (color scale) with significance contours from H.E.S.S. observations of W28 
overlaid (green). The black contours follow the 95 % and 99 % confidence level re- 
gions for the unidentified EGRET source GRO J1601-2320 ( |Aharonian et al.|2008e| l. 
Right panel (b): The color scale indicates the smoothed excess map obtained from 
MAGIC observations of IC 443. The molecular gas density (CO) is traced with the 
cyan colored contours, X-ray contours (purple) from ROSAT observations, gamma- 
ray contours (EGRET) in black. See captions of Fig. 1 in Albert et al. (2007c]) for 
references. 



Galactic SNR are of MM supports this scenario. 
The observation of VHE gamma-rays from MM SNR like W28 (|Aharonian et al. 



|2008e| l and IC443 ( |Albert et al.|2007c| is - independent of the underlying radiation 
mechanism - direct evidence for particle acceleration to multi-TeV energies in these 
objects (see Fig. 13}. However, one would expect that the evolved and slow shock 
present in these objects is not strong enough to accelerate particles at the current 
stage of evolution. Furthermore, the observed Gamma-ray emission seems to coin- 
cide well with regions of dense molecular gas. In the north-eastern region of W28, 
some of the clouds must have already been interacting with the expanding blast wave. 
The interaction can be traced by maser emission in the shocked gas (see Fig. [3^). The 
observation of VHE-emission from IC 443, another MM type SNR, supports a sim- 
ilar scenario of an evolved SNR interacting with dense and cold molecular gas (see 
Fig. J3]}). In this scenario, electrons can not produce the observed VHE-emission 
because of the rapid cooling through inverse Compton and synchrotron emission, as 
well as Bremsstrahlung in the dense gas. Typical cooling times of electrons at multi- 
TeV energies would be of the order of a few hundred years and by far too short to 
allow for a sizeable population of electrons to be present in an evolved SNR with- 
out ongoing acceleration. An alternative and much more favorable interpretation 
for MM SNRs requires that accelerated nuclei have been partially confined to the 
SNR and its environment and produce efficiently gamma-rays through 7T°-decay in 
the dense target material of the molecular clouds. Gabici & Aharonian (2007) have 
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investigated the effect of cosmic rays that have left the accelerator and produce an 
observable VHE emission in a nearby cloud of molecular gas. 



3.3 Pulsar wind nebulae and composite supernova remnants 

While shell type SNR are considered to be the likely source of the Galactic cos- 
mic rays, the blind survey of the Galactic plane (Aharonian et al. 2005b 2006k) 
performed with the H.E.S.S. telescopes has revealed a comparably small number of 
shell-type SNR while VHE emission from pulsar wind nebula (PWN) systems has 
been observed more frequently (Gaens kTr & Slane|2006 l. In total, 20 VHE sources 



(including candidate associations) have been discovered which are very likely pow- 
ered by isolated pulsars. 

The population of "TeV-Plerions" ( |Horns et al.||20061> includes young objects lik e 



Kes-75/PSR J1846-0258 (spin-down age X = 723 yrs) ( | A. D jannati- Atai et al.|2007) , 



MSH 15-52 ( |Aharonian et al.|2005c) driven by PSR J15 13-5908 (t = 1.55 kyrs) as 



well as evolved systems like the PWN driven by PSR B1823-13 (T = 21.4 kyrs) as- 
sociated with HESS J1825-137 ( |Aharonian et al.|2005j] ). Apparently, PWN systems 
are active VHE sources for a few 10000 years thus exceeding the time that a SNR 
evolves through the Sedov-stage of shell type SNR during which presumably VHE 
emission is most likely produced. The longer time during which PWNe are active 
particle accelerators explains the larger number of PWN systems active at any given 
time in the Galaxy in comparison to shell-type SNRs. 

However, the contribution of PWN to the Galactic cosmic rays remains an open ques- 
tion as long as the nature of the particles accelerated in these systems is not clarified. 
Currently, for most of the VHE emitting PWNe, a leptonic origin of the VHE emis- 
sion can qualitatively explain the multiwavelength morphology and energy spectra. 
As an example, HESS J1825-137 has been studied intensively both in X-rays as well 
as in VHE-gamma-rays where the source extension appears to be larger by a factor 
of w 6 than in the X-ray band. This larger size can be explained by a greater loss 
time of electrons radiating VHE gamma-rays via inverse Compton scattering than 
the loss time of electrons radiating synchrotron X-rays (|Aharonian et al. 2006h). 



The energetic electrons injected by the pulsar PSR B 1823- 13 loose rapidly energy 
in the high magnetic field environment close to the pulsar where the X-ray nebula 
is observed (Gaensl er et al.|2003) . While the particles leave the high magnetic field 
environment, they presumably radiate only at wave-lengths longer than the X-ray 
band in the synchrotron channel while inverse Compton scattering off the cosmic 
microwave-background is the dominant energy loss process leading to observable, 
extended VHE emission. In this model, the slowly decaying spin-down power of the 
pulsar is accumulated in an expanding bubble containing the TeV-emitting electrons. 
This scenario is supported by softening of the observed VHE spectra with increasing 
distance to the compact X-ray PWN system (Aharoni an et al.|20 06h). 
The observed morphology of the PWN is often asymmetric and not centered on the 
pulsar. While in some cases, the velocity of the pulsar can accommodate for the off- 
set, other objects require an alternative explanation (e.g. Vela X): here, the asymmet- 
ric reverse shock of the SNR shell interacts with the relativistic PWN gas "crushing" 
and pushing it off-centered ( |Blondin et al.|200T Ferreira & de Jager 2008 ). 
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A contribution to the observed VHE gamma-ray flux from nuclei accelerated by the 
pulsar is not excluded but seems to be less important than the contribution of elec- 
trons for most systems. 

The Vela-X PWN system may be an exception. Here, the spin down power of the pul- 
sar is by far larger than the acceleration rate of electrons required to drive the PWN 
system. In this sense, there is missing energy that can be naturally explained by a 
nucleonic component in the wind that drives the acceleration of electrons through 
resonant wave absorption ( Arons & Tavani|1994[ ) and produces the bulk of the ob- 
served VHE emission ( Aharonian et al.|2006i| l. The VHE spectrum can be explained 
quite well in terms of over-all energetics (thus solving the problem of "missing" en- 
ergy) as well as its shape ( |Horns et al.||2006) . A crucial test will be the detection 
of high energy neutrinos from this object that may be the brightest steady neutrino 
source in the sky (Kap peset al.|2007b|a I. 

Finally, composite systems where a plerionic, non-thermal X-ray component is ob- 
served to be embedded in a radio-emitting shell like GO. 9+0. 1 have been found to 
emit VHE gamma-rays (Aharonian et al. 2005g). In addition to known composite 
SNRs, VHE observations have helped to discover and identify new systems: X-ray 
observations of HESS J1813-178 ( |Aharonian et al.|2005b) > carried out with XMM- 
Newton have revealed a previously unknown extended X-ray source inside the radio- 
shell of SNR G12. 82-0.02 ( |Funk et al.||2007] > putting this object in the category of 
composite SNR systems. Most likely, the VHE-emission is produced in the central 
PWN systems while the radio-shell is VHE-dim. However, it should be noted, that 
the spatial extension of these two systems is barely resolved at VHE energies. 



3.4 X-ray binary systems 

Before the imaging air Cherenkov technique was widely established, several groups 
had claimed a number of X-ray binary sources to be periodic (e.g Her X-l, Cyg X- 
3). Cyg X-3 had been claimed to be a transient source of energetic particles up to 
PeV energies ( Samorski & Stamm|1983 i, for a review of early claims of detections 



see |Weekes ( 1992} . With the notable exception of a possible burst observed from 



the direction of GRS 1915+105 ( |Aharonian & Heinzelmann|1998| l, none of the XRB 
systems had been confirmed until recently to be VHE gamma-ray emitters. Finally 
4 XRB systems listed in Table|2]have been detected to be VHE gamma-ray sources. 
All of these objects show variability in the VHE band either linked to the orbital 
motion as for LS I +61 303 (|Albert et al.|2006f]paierp007ll, LS 5039 ( |Aharonian| 



[eTaLl|2006c| l, and PSR B1259-63 (Aharonian et al. 2005d), or transient activity as 
for Cyg X- 1 ( [Albert et al.|2007h) . 

Given that these objects are currently the only variable Galactic sources of VHE 
gamma-rays, the interpretation of their properties is quite different from the other 
source types. The orbital parameters as well as the physical conditions in these sys- 
tems are reasonably well-known and allow for a detailed modelling of the processes 
leading to acceleration and VHE emission. In this sense, one can consider the sys- 
tems as "laboratories". 

The pulsar PSR B 1259-63 (P = 48 ms) is in a highly eccentric orbit (f rbit = 
3.4 yrs) around the Be-type companion star SS2883. VHE emission is predomi- 
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Table 2: X-ray binary systems detected with imaging air Cherenkov telescopes.. 



1UC11L111C1 


Van ahi 1i tv 


rj 
[kpc] 


If 

[10 33 ergs/s] 


oys lciii 


PSRB 1259-63 


P=3.4a 


1.5 


1.6 


PSR(47 ms)/B2e 


LS 5039 


P=3.9d 


2.5 


8.7 


BH or NS/06.5V 


LS I +61 303 


P = 26.5d 


2 


2.5 


NS/BOVe 


CygX-1 


r var ~ 80 min 


2.2 


1.6 


BH/09.7 lab 



"Integrated from 1 to 10 TeV 

*BH: Black hole candidate, NS: Neutron star, PSR: Pulsar 



nantly produced around the periastron where the effects of adiabatic and radiative 



cooling due to inverse Compton as well as synchrotron emission are strongest ( Aha- 



|ronian et al.|2005d I. The IC component had been predicted successfully prior to the 
detection ( |Kirk et al.|1999| l, albeit the observed temporal behavior of the source as it 
passed through the periastron in 2004 was quite different from the expectation ( |Aha- 



ronian et al.|2005d ). Refined models including a more detailed treatment of inverse 



Compton scattering in energetic and anisotropic radiation fields as well as adiabatic 
energy losses can reproduce the observed light curve and provide some crucial pre- 
dictions to be probed with new observations (Khang ulyan et al.|2007[ ). Alternative 
models invoking the presence of a nucleonic component that leads to VHE emission 
through the interaction with the disk-outflow of the Be-star have been proposed as 
well ( |Kawachi et al.|2004| |Neronov & Chernyakova|2007 1. 

Unfortunately, the periastron passage in 2007 took place during the time when PSR 
B 1259+63 culminated during the day time. A crucial observation will be only possi- 
ble during the periastron passage in 201 1 when the system will be visible for ground 
based VHE instruments as well as for Fermi. 

The two systems LS 5039 as well as LS 1 +61 303 have shorter orbital periods 
of 3.9 and 26.5 days, respectively. Initially, both objects have been considered to be 
micro-quasars with a steady radio-jet feature. There is an ongoing debate whether 
this picture may have to be modified. Radio observations of LS 1+61 303 indicate 
that the orientation of the "jet" changes during the orbit (Dhawa n et al.|20 06 ) which 
is more in line with the cometary tail predicted by Dubus (2006b ). The cometary 
tail is a consequence of the interaction of a pulsar wind with the stellar outflow. 
For LS 5039, the current high resolution radio observations also indicate a change 
in the orientation of the jet-like structure (Ribo et al. 2008). However, since the 
observations do not cover a full orbital period, the results are not conclusive yet. 
Given the uncertainty of the nature of the compact object, the VHE emission from 
LS 5039 has been interpreted in the context of a micro-quasar model (Bosch-Ramon 
et al.|2006| ) as well as in an interacting PWN model (Dubus 2006b]). 



A crucial consideration is the absorption of VHE photons due to pair-production 
in the photon field of the stellar companion. If the line-of-sight towards the VHE 
emission region passes within roughly 10 14 cm of the companion star of LS 5039, 
absorption will produce visible effects (Dubus 2006a) varying with the phase of the 
orbit. The absorption effect is energy dependent and should lead to a pronounced 
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E(eV) 



Figure 4: The VHE energy spectrum of LS 5039 for two distinct orbital phases. The 
INFC corresponds to the inferior conjunction, while SUPC is the phase interval at 
the superior conjunction with the periastron passage ( Aha ronian et al.|2006c[ >. 



variability at energies between 200-400 GeV. However, the observations indicate 
that the flux at 200-400 GeV remains almost unchanged along the orbit while the 
modulation is most pronounced at TeV-energies (see Fig. |4jl. 

Secondary particle production in cascades can in principle reproduce the observed 
variability (Bednarek 2007b) as well as models where the emission takes place at a 
larger distance to the stellar companion (Khangulyan et al. 2008 ). 



3.5 Star-forming regions 

So far, the objects under scrutiny are related to the final stages of stellar evolution 
and their remnants. However, gamma-rays are also produced in molecular clouds 
which are the cradles of star formation and possibly during the final mega-year of 
stellar evolution where massive stars (> 15 M Q ) drive powerful winds. The wind- 
phase of stellar evolution is characterized by mass loss rates of up to 10~ 5 M /year 
and terminal velocities of a few 1000 km/s. These early-type stars are usually born 
in open stellar associations in the spiral arms of the Galactic disk. The massive O 
stars in these associations often evolve into the Wolf-Rayet phase where strong stel- 
lar winds produce a hot tenuous cavity in the interstellar medium. In such systems 
up to 10 39 ergs/s are dissipated in the form of kinetic energy from stellar winds, even 
before member stars evolve into supernovae. Provided that shocks will form either 



in the termination region of cumulative winds (Domingo-Santamana & Torres 2006) 
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or in wind- wind interaction regions ( Reimer et al.|2006 De Becker|2007 1, some of 
this power can be converted into cosmic-ray acceleration. 

The detection of VHE-emission from Westerlund-2 is a crucial step towards estab- 
lishing the role of stellar driven cosmic ray acceleration (see Fig . [5] from | Aharonian 



|et al.| ( [2007d| l). Westerlund-2 is a young star forming region with roughly 20 O-type 
stars and the Wolf-Rayet systems WR 20a&b ( |Rauw et al.||2007[ ). The total power 
released in the stellar winds is estimated to be 5 x 10 SI ergs/s. The observed VHE 
emission can be easily produced if roughly 1 % of the kinetic energy is converted 
into the acceleration of electrons (Manolakou et al.|2007)|Reimer et al.|2007] ). 
Prior to the discovery of VHE emission from Westerlund-2, the first unidentified 
VHE source TeV J2032+4130 had been associated with the remarkably powerful 
OB association in the Cygnus arm. The HEGRA discovery of TeV J2032+4130 
(Aharoni an et al.|2002al |2005a|> has b een recently confirmed by MAGIC ( |Albert & 
|for the MAGIC Collaboration||2008|l as w ell as by the MILAGRO and VERITAS 
groups ( |Abdo et al.||2007a|b| |Lang eTaLl[2004) [Konopelko et al.|[2007] l. Remark- 
ably, observations with the MILAGRO detector show an extended (« 3°) source of 
> 10 TeV photons centered on TeV J2032+4130 ( |Abdo et al.|2007a) . This source has 
received considerable multi-wavelength coverage including deep radio observations 
which show indications for extended radio emission coinciding with the VHE source 
( Paredes et al.|20 07 ). Recent X-ray observations with XMM-Newton have revealed 
faint non-thermal extended X-ray emission co-located with the VHE-source ( |Horns| 
|et al.|2007) . 

Future observations of this region with the VERITAS telescope array will be very 
helpful to disentangle possible multiple-sources and to study the source size at dif- 
ferent energies. 

Besides Westerlund-2 and Cyg OB2, the MILAGRO source MGRO J20 19+37 ( |Abdo| 
|et al.|2007a[ l has been associated with the young stellar association Berk-87 includ- 
ing the Wolf-Rayet object WR 142 ( |Bednarek|2007a| l. 

A systematic search for VHE emission with the HEGRA telescopes from young open 
clusters has so far produced only upper limits (Aharon ian et al.|2 006a). 



3.6 Unidentified Galactic sources 

The sources discussed in the previous section have been identified or are associated 
with known objects. However, a number of objects remain to be identified. A total 
of 10-20 objects have no clear or even multiple counterparts. Among these objects, 
the Galactic center source and the point-like source in the Monoceros region are 
highlighted in the following. For a recent summary of unidentified sources, see e.g. 
[Funk| ( |2007l l; |Aharonian et"aT] | |2008c) . 



Galactic center region. The central r egion of our Galaxy contains a super-m assive 



black hole with a mass of 3 x 10 6 M JSchodel et al 



2002 



Ghez et al. 



2005 ). Given 



the faintness of this object, the accretion rate must be < 10 of the Eddington-rate. 
The mm-, IR-, and X-ray sources associated with Sgr A* have been observed to show 
variability including flares and outbursts with a time-scale of hours ( |Genzel et al.| 
2003 ). The shortness of the flares constrains the emission region to be compact with 
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Right ascension (J2000) 



Figure 5: Westerlund 2 (Aharonian 
|et al.||2007d| ): The color scale indi- 
cates the smoothed VHE excess map 
(see the bottom left inlaid figure for 
the point-spread function), the con- 
tours are the significance levels start- 
ing at 5 <T and increasing by 2. The 
filled circle and error bars mark the 
position and uncertainty of the VHE 
source (1 a). The dashed circle marks 
the position and extension of Wester- 
lund 2, the upright triangle is at the 
position of the Wolf-Rayet binary sys- 
tem WR 20a and the downward trian- 
gle is at the position of WR 20. 



a spatial extension r < f var -c 10 14 cm (t var /hr) « 100r G (M/3 • lO 6 M ) _1 f V£ , r /hr. 
Observations with the CANGAROO III telescopes as well as at large zenith angles 
with the Whipple 10m telescope revealed a VHE gamma-ray source co-located with 
the Galactic center (Tsuchiya et al. 2004, Kosack et al. 2004). The positional ac- 
curacy as well as spectral measurements have been considerably improved with the 
observations of the H.E.S.S. -telescopes ( | Aharonian et al.|2004b") . The MAGIC tele- 



scope has confirmed the earlier findings in independent observations (Albert et al. 
|2006dD . 

The detection of a point-like source with a hard power-law spectrum was shortly 
after the discovery discussed in the context of emission scenarios near to the super- 
massive black hole (at distances larger than w 10 r<j in order to avoid internal ab- 
sorption) from e.g. ultra-high energy protons emitting synchrotron and curvature 
radiat ion in the high mag netic field s that are believed to be threaded in th e accretion 
disk ( Aharonia n & Neronov|2005j l or 7T°-decay from energetic protons (Ballantyne 



|et al.|200"7{|Liu et al.|2006|l. A t a larger distance from Sgr A*, acceleration in SNRs 
(Erly kin & Worfendale|2007[ l including Sgr A E ast (|Crocker et al.|2 005 ) as well as 
acceleration in stellar winds have been proposed (Qua taert & Loeb|2005fr . 
The initial data were consistent with a Dark Matter annihilation scenario albeit re- 
quiring uncomfortably large masses beyond 20 TeV for the annihilating particles 
( |Horns|2005| l. 

With more data, the H.E.S.S. experiment has constrained the energy spectrum to ex- 



tend beyond 10 TeV and to deviate from a WIMP-annihilation scenario (Aharonian 
2006j i 14 . With the increased exposure, faint spatially extended VHE emission 



et al. 



along the Galactic ridge was detected (Aha ronian et al.|2"00 6g). The accuracy of re 



construction the position of the VHE-emitting point source has been improved to the 

14 see howevei jBringmann et al.|{2008) for a discussion on the effects of internal Bremsstrahlung pho- 
tons which are important for heavy WIMP annihilation in order to reduce the helicity suppression for the 
annihilation channel 
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Figure 6: Left panel: Chandra X-ray and SINFONI near-IR composite from | Wang] 
|et al.| ( (2006l > with additional annotations. The circular region corresponds to the com- 
bined systematic and statistical uncertainty for the HESS-source associated with the 
Galactic center. Right panel: X-ray and VHE-gamma-ray light curve from simul- 
taneous observations of Sgr A*. The H.E.S.S. light curve has been scaled to match 
the Chandra light curve in units of counts per second. Note the comparably small 
error bars of the H.E.S.S. observations - a flare of similar magnitude as observed in 
the X-ray band would have been clearly detected. 



level of the systematic pointing uncertainty of ss 6 arc seconds ( |van Eldik et al. 2007 1 



and excludes Sgr A East. The error box for the VHE source still encompasses at least 
three possible candidates for VHE emission (see Fig. [6^): Sgr A*, a low-mass X-ray 
binary system (LMXRB) ( |Muno et aT1|2005| l, the stellar system IRS 13 ( |Baganoff| 
|et al.|2003) , and the PWN G359.95-0.04 ( |Wang et al.|2"006| . The PWN system has 



been argued to be a possible candidate to explain the observed VHE emission (Hin 
|ton & Aharonian|2007| ). 

An association with a variable source like Sgr A* would be supported if simulta- 
neous multi-wavelength (MWL) observations of Sgr A* during an outburst would 
show a correlation between e.g. X-ray and gamma-ray flux. During a dedicated 
MWL observation of Sgr A* in summer 2005 with H.E.S.S. and the Chandra X-ray 



telescope, a typical X-ray flare was detected with the Chandra instrument dHmtan 



et al.|2007a i. The simultaneously taken H.E.S.S. -data show however no variability 
even though the observed rate of VHE photons is sufficiently large to detect a flare 
of similar strength as seen in the Chandra data (see Fig.|6|)). The absence of corre- 
lation would disfavor at least a common origin of the X-ray and VHE emission from 
Sgr A*. 

An association of the VHE source with the LMXRB could in principle be demon- 
strated by correlated activity in different wavelength-bands. However, so far no such 
observations are available. 

At this point it is not possible to discern between the different source candidates 
nor exclude a dark matter origin. 
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Figure 7: Monoceros source ( A haronian] 
|et al.||2007h| ): The grey-scale illustrates 
the velocity-integrated (0-30 km s _1 ) 
12 CO emission as measured with the 



NANTEN telescope ( |Mizuno & Fukui 
2004). The yellow-contours trace the 
4 and 6 <7 level of detection with the 
H.E.S.S. telescopes. The cyan contours 
indicate the 8.35 GHz radio measure- 
ments from Langston et al. (2000) while 
the green contours mark the 95 % and 
99 % confidence regions for the position 
of the EGRET source 3EG 0634+0521. 
SAX J0635. 2+0533 is a pulsar binary sys- 
tem. The positions of Be-stars are indi- 
cated with pink stars. 



Monoceros source. The newly discovered population of unidentified VHE sources 
in the Galactic plane form a quite homogeneous sample of objects with respect to 
their spatial extension and energy spectra. It should be emphasized however, that the 
similarity of the VHE characteristics does not necessarily imply that the underlying 
sources are similar. Certainly, observational biases are introduced by e.g. the limited 
sensitivity for detecting sources with low surface brightness. 

Among the unidentified sources, there is - besides the Galactic center - only one 
more unresolved object: HESS J0632+058 which is close to the rim of the Mono- 
ceros Loop SNR (Aharonian et al. 2007h). The radial extension of this object is 
constrained to be smaller than 2' (95% c.l.). No strong indications for variability 
are observed. The energy spectrum of the source is well-described by a power-law 
(dN/dE °< £-2.5±0.3) ^yjjjQjj j s a t so fj- enc [ f the distribution of energy spectra 
observed from unidentified VHE-sources. 

For an unresolved source (r < 2'), the search for a counterpart is simpler than for ex- 
tended objects. In Fig. [7] the position of HESS J0632+058 is combined with multi- 
wavelength observations. Only a few candidate objects coincide spatially with the 



VHE source including an unidentified EGRET source 3EG J0634+0521 (Hartman 



et al 



1999) ), an unidentified X-ray source 1RXS J063258. 3+054857 ( |Voges et al 



2000), and a BOpe type star MWC 148 which may be a binary system similar to 
PSR B 1259-63. Alternatively, it could be a new type of VHE source associated with 
an isolated stellar system. 

3.7 Extra-galactic sources 

The blind survey of the Galactic plane has revealed an interesting and previously 
unknown population of VHE-sources. Unfortunately, a survey with similar sensitiv- 
ity of the extra-galactic sky has to be postponed until wide-field-of-view Cherenkov 
telescopes or the next generation of Water-Cherenkov-experiments like HAWC will 
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come online (see also next section). Currently, the observations have focussed on 
Blazars ("BL Lac type Quasars"). 



Historically, the discovery of VHE emission from the Blazars Mkn 421 (Punch et al. 



1992) and Mkn 501 (Qui nn et aL[T9 96 ) have opened the field of TeV-Blazar obser- 



vations. These objects are in general characterized by a featureless optical spectrum, 
non-thermal X-ray emission, and violent variability in the optical wavelength-band 
as well as in X-rays. The spectral energy distribution is dominated by two broad max- 
ima - one in the optic al-to-X-rays and one in the gamma-ray-to- VHE-gamma-ray do- 
main. The low-energy maximum is usually assumed to be produced by synchrotron 
emitting electrons while the high energy bump would be due to inverse-Compton 
scattering. 

TeV-Blazars can be considered "clean" systems where the main source of seed pho- 
tons are actually the synchrotron photons emitted by the electrons themselves, com- 
monly called Synchrotron-Self-Compton (SSC) mechanism. Within the unified scheme, 
Blazars are believed to be Fanaroff-Riley Type I galaxies with a jet axis pointing 
close to the line-of-sight. 

In a simple one-zone-SSC model, the VHE luminosity is closely linked to the syn- 
chrotron luminosity which in turn is related to the particle density and provides 
seed photons for inverse Compton scattering. A simple but very successful selec- 
tion scheme of possible VHE emitting candidates is based on the X-ray or radio flux 
of Blazars (Costamante & Ghisellini 2002). Most of the Blazars that have been pre- 
dicted to be detectable according to this scheme have been actually discovered by 
the H.E.S.S. and MAGIC collaboration. 

With the growing number of known Blazars at different red-shifts (see Table [3] for 
a list of extra-galactic VHE-gamma-ray sources), a number of issues and questions 
can be addressed: 

• Is a simple one-zone SSC model sufficient to explain all (multi-wavelength) 
data including correlated variability patterns or is a more complicated model 
required (stratified jets, multiple zones model)? 

• Where does the acceleration take place? 

• What is the underlying acceleration mechanism in Blazars and what is the 
composition of the jet? This is a fundamental question related to the jet physics 
and the energetics of AGN. 

• What is the duty cycle of Blazars? 

• How strong is the beaming? In the relativistic outflow, a co-moving isotropic 
emission is strongly beamed in the forward direction. With increasing bulk 
Lorentz factor T, the opening angle gets narrower <x r _1 . The higher the 
beaming, the smaller the likelihood that we can observe the emission. In con- 
sequence, a high beaming factor can lead to dramatically increased number of 
Blazars that are pointing not in our direction. 

• Is there internal absorption (due to photon-photon scattering with pair-production)? 

• How strong is the absorption on the extra-galactic background light (EBL)? 
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In the following, I will highlight observations from Blazars as well as from non- 
Blazar VHE sources including M87 and the flat-spectrum-radio quasar 3C279. 



3.8 Blazars 

A very spectacular feature of VHE emission from Blazars is the short time-scale 
variability. The most impressive examples of short flares have been observed in 
2006 from the nearby object PKS 2155-305 (red-shift z = 0.116). This object was 
the first Blazar discovered in the southern hemisphere by the Durham group using 
their Mark VI telescope in Australia (Chadw ick et al.| 19 99 ). The H.E.S.S. telescopes 
observed their "first gamma-light" from this source already in 2002 when the first of 
the four H.E.S.S. telescopes went online ( [Punch ||2007] l. Since PKS 2155-305 is a 
remarkably bright object, it can be easily detected in quiescent state with the system 
of four H.E.S.S. telescopes. Detailed measurements of the energy spectrum show 
despite widely different flux states only marginal deviations from a soft power-law 
with photon index sa 3.3 ( Aharonian et al. 2005f ). PKS 2155-305 has been the target 
of a number of simultaneous multi-wavelength-campaigns together with the RXTE 
X-ray satellite as well as in conjunction with the Chandra, XMM-Newton, and the 



SWIFT X-ray telescopes ( Aharonian et al.|20 05i). During summer 2006, the source 



went into an unprecedented high state reaching flux levels more than a factor of 
100 higher than the quiescent flux (Ah aronian et al.|20 07c). During the peak flares, 
the rate of detected VHE gamma-rays reached values exceeding 1 Hz which marks 
the highest rate ever detected at VHE energies. During the night with the highest 
flux, a sequence of flares can be discerned with rise and decay timescales of minutes 
(see Fig. [HJ. This immediately constrains the emission region to be very compact 
and to be moving with a large bulk Lorentz factor to overcome pair-opacity and 
to relax the size constraint (Begelman et al. 2008). Under the assumption that the 
variable emission is produced close to the super-massive black hole, two interesting 
conclusions have been pointed out recently: (a) The variability time scale constrains 
the mass of the black hole to be of the order of 10 7 M® roughly 2 orders of magnitude 
smaller than the general correlation of the mass with luminosity of the host galaxy 
suggests ( |Celotti et al.][T99"8| |Neronov et aTJ|2008l |Bettoni et al."1|2003j >. (b) The 
apparent periodicity of the subflares observed could potentially provide a measure of 
the angular momentum of the black hole ( Neronov et al.|20 08 ). 



3.9 Other active galactic nuclei 

Besides the Blazars, only two non-Blazar extra-galactic objects have been discovered 
so far: M 87 and 3C 279 15 



M87 Besides the collective study of a Blazar sample (Wagner 2008), it is interest 



ing to closely examine nearby FR-I radio galaxies as these objects are considered to 

15 A signal from Centaurus A at energies above 300 GeV had been claimed on the basis of non-imaging 
observations of air Cherenkov light at the level of 4.5 o~ i Grindlay et al. 1975 1. This claim awaits confir- 
mation. 
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Table 3: List of extragalactic VHE-sources sorted by red-shift. Whenever available, 
different source activity states are characterized . 



Source 


Type 


z 


d L " 


l°gl0 W 


^measured 


Ref. 


M87 


FRI 


0.0044 


23 


41.1 


2.6±0.4 


C 










41.6 


2.2±0.2 


C 


Mkn421 
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0.030 


130 


44.8 


3.0±0.2 


a 










45.3 


2.06 ±0.03 


d 


Mkn501 


HBL 
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142 


44.3 


2.45 ±0.07 


e 










45.2 


2.09 ±0.03 


e 


1ES 2344+514 


HBL 


0.044 


183 


43.9 


2.9±0.2 


f 










45.2 


2.5±0.2 


8 


Mkn 180 


HBL 


0.045 


194 


44.0 


3.3±0.7 


h 


1ES 1959+650 


HBL 


0.047 


198 


44.2 


2.7±0.1 


i 










44.9 


1.8±0.2 


j 


BL Lacertte 


LBL 


0.069 


293 


44.0 


3.6±0.5 


k 


PKS 0548-322 


HBL 


0.069 


300 


43.6 


2.8±0.3 


I 


PKS 2005-489 


HBL 


0.071 


306 


44.2 


4.0 ±0.4 


ni 


RGB J0152+017 


HBL 


0.080 


345 


44.0 


3.0±0.4 


a 


W Comae 


IBL 


0.102 


452 


44.9 


3.8±0.4 


N 


PKS 2155-304 


HBL 


0.116 


515 


44.1 


3.4±0.1 













46.8 


2.7±0.1 


P 








Break at 0.43 TeV 


3.5 ±0.1 


P 


H 1426+428 


HBL 


0.129 
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45.9 


3.7 ±0.4 


q 


1ES 0806+524 


HBL 


0.138 
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43.8 


assumed 2.6 


Q 


1ES 0229+200 


HBL 


0.139 


633 


44.3 


2.5 ±0.2 


r 


H 2356-309 


HBL 


0.165 
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44.5 


3.1 ±0.2 


s 


1ES 1218+304 


HBL 


0.182 
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45.2 


3.0 ±0.4 


t 


1ES 0347-121 


HBL 
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44.8 


3.1 ±0.3 


a 


1ES 1101-232 


HBL 
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877 


44.8 


2.9±0.2 


V 


1ES 1011+496 


HBL 


0.212 
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45.5 


4.0±0.5 


w 


PG 1553+113 


HBL 


>0.3 


> 1498 


>46.1 


4.5 ±0.3 


X 






<0.74 


<4437 


<47.1 




X 


S5 0716+714 


IBL 


0.31 


1557 


46.1 


assumed 3.8 


XY 


3C279 


FSRQ 


0.536 


2998 


46.7 


4.1 ±0.7 


y 



0.73, H Q = 73km s^'Mptr 1 

1 integrated betwee n 0.1 and 1 TeV 

A lbert et al j|2007g[, ^Aibert et al.](2 "007e| 
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Time - MJD53944.0 [min] 

Figure 8: The light curve from PKS 2155-304 in the night July-28 2006 as measured 
with H.E.S.S. above 200 GeV. The data are binned in 1-minute intervals. For com- 
parison, the horizontal dotted line represents I(>200 GeV) observed from the Crab 
Nebula. The curve is a fit of a superposition of five individual bursts and a constant 
flux. See |Aharonian et al. (2007c) for additional details. 



be mis-aligned Blazars. The closest FR-I objects are M87 and Cen-A. VHE gamma- 



ray emission from M87 had been initially reported from the HEGRA group ( Aharo- 



nian et al. 2003b ) and later verified by the Whipple collaboration (Le Bohec et al. 
2004| ), as "well as by the H.E.S.S. group ( Aharonian et al.|2006b ) and most recently 



with the VERITAS telescopes ( |Acciari et al.|2008b| ). 
While the initial sensitivity was only sufficient to detect a signal and coarsely con- 
strain the energy spectrum, following observations mainly with the H.E.S.S. tele- 
scopes have lead to a number of important discoveries: the energy spectrum of M87 
extends to energies beyond 10 TeV and even more spectacular, the observed emission 
is highly variable on time-scales as short as days (Aharonian et al. 2006b). The hard 
spectrum is a surprise and disfavors the models suggested by |Georganopoulos et al 



(2005) and |Reimer et al.| ( |2004[ ) where the predicted gamma-ray spectrum is soft as 
a consequence of the assumption of a small Doppler-factor derived from radio mea- 
surements of the jet dynamics. It should be noted however, that also other nearby 
Blazars that have been observed with high resolution radio interferometers do not 
show strong super-luminal motion indicative for large Doppler-factors. It seems that 
there is growing evidence for a discrepancy between high Doppler factors implied 
by modelling the VHE emission of Blazars with single zone SSC model and the lack 
of evidence for strong super-luminal motion in these objects from radio observations 
(Piner et al. 2008[ ). The fast variability observed from M87 leads to interesting impli- 
cations. First of all, the variability excludes models where the VHE gamma-rays are 
produced in processes like Cosmic ray interactions in the host galaxy (Pfrom mer &| 
EnBlin 2003 ) or even Dark Matter annihilation (Baltz et a l.|2000[ ). These processes 
should not lead to temporal variation of the observed emission. Given that the vari- 
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Figure 9: From Acciari et al. ( 2008b J: The left panel combines the VHE and X- 
ray light curves from observations of M87 with ground based Cherenkov telescopes 
as well as with satellite instruments on-board the RXTE satellite (ASM) and the 
Chandra observatory. The right panel correlates the yearly VHE fluxes with the 
appropriately averaged ASM fluxes. The dashed curve is a linear correlation and the 
solid curve a squared function. 



ability time scale can not be smaller than the light crossing time, the constraint on 
the size of the emission region is severe and excludes regions in the outer radio jet as 
well. A possible candidate for the emission could be the notoriously variable HST-1 
knot located less than 1 arc sec off the core. The most recent long-term light-curve 
obtained with the all-sky-monitor on-board the RXTE satellite as well as a dedicated 
Chandra monitoring seems to support a correlation of the HST-1 activity in X-rays 
and the observed VHE emission (see Fig. [9]). 



3C279 This object belongs to the so-called flat-spectrum radio quasars (FSRQs). 
These objects are long-known to be emitters of GeV gamma-rays that have been de- 
tected with EGRET (Har tman et al. 1 1 992] |200 1 ) > . There are two aspects which single 
out this source: (a) 3C279 has a comparably large red-shift of z = 0.536 and (b) it is 
a new type of AGN which shows acceleration to TeV energies. 
The EGRET Blazars are mostly members of the FSRQ-class and are characterized 
by a spectral energy distribution with a peak in the EGRET energy range, dominat- 
ing the total luminosity. The observed emission is commonly attributed to external 
Comptonization of photons from the broad-line region and possibly from the accre- 
tion disk. Given that FSRQs have intrinsically a higher gamma-ray luminosity than 
Blazars, they can be in principle observed to higher red-shifts opening up the possi- 
bility to finally detect absorption features in the VHE gamma-ray spectra (see also 
next section) of extra-galactic objects (see also next section). The detection of VHE 
emission from 3C279 with the MAGIC telescope (Teshima et al. 2007[ ) indicates 
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variability which seems to be correlated with changes in the optical flux. Further 
observations of this source and a measurement of the VHE energy spectrum could 
provide important clues on the level of extra-galactic background light (see next sec- 
tion). Even in the case of a low level of EBL, only 10~ 3 of the emitted flux will be 
observable at 1 TeV ( |Raue & Mazin|2008| l. 



4 "Secondary physics": Propagation of gamma-rays 
4.1 EBL absorption 

VHE emission from AGN appears to be more common than what could have been 
hoped for and the temporal variability observed at the highest energy is providing us 
an important and unique in-sight into the "engine-room" of AGN. The understand- 
ing of processes leading to VHE-emission hinges however on a good understanding 
of pair-absorption processes of the type YvheYebl — > e + e~ . This process has been 
known for a long time to be relevant for the propagation of TeV-photons over cos- 
mological distances ( [Gould & Schre der 1966). 



With the growing number of Blazars observed at various red-shifts it has become 
feasible to use VHE -photons as probes to constrain or even measure the background 
photon density in extra-galactic space. Until now, it has however not been possible 
to clearly establish absorption in the measured Blazar spectra because the intrinsic 
emission spectrum is not very well known and may suffer additional internal ab- 
sorption (Aharonian et al.|2008d| ). When invoking a constraint on the shape of the 
intrinsic spectrum, it is however possible to derive (model-dependent) upper limits 
on the EBL density ( |Dwek & Krennrich||20"05] |Aharonian et al.||2006el |Mazin &| 
|Raue|2007) . The situation will be greatly improved when the first broad-band spec- 
tra of Blazars will be obtained that should show a transition from the optically thin 
to the optically thick part. The detection of such an absorption feature will in turn 
lead to a detection of the EBL. In combination with increasing statistics and cover- 
age over varying red-shift, it will be feasible to investigate the evolution of the EBL 
and finally, derive in a completely independent way constraints for cosmological 
parameters (Blanch & Martinez 2005a b). 



4.2 Lorentz invariance violation 

VHE photons are useful probes for Lorentz invariance violating (LIV) processes. 
A number of theoretical approaches towards a quantum theory of Gravity predict a 
structure of space-time at the Planck scale. This will result in a modified disper- 
sion relation for photons of the form c 2 p 2 — E 2 (l + 1 (E/Mqq) + 0(E /Mq G ) 2 ) see 
e.g. [Sarkar (2002]) for a review. This dispersion relation affects the time-of- flight of 
photons of different energies and the kinematics of interaction, suppressing e.g. the 
pair-production process of photons. 

Currently, fast flares from Mkn 421 ( |Gaidos et aL] [T996 ), Mkn 501 ( |Albert et aT] 



2007g), and PKS 2155-304 (Ahar onian et al.|2007c| ) have been observed. The latter 



is potentially the most constraining observation as the red-shift is larger and the vari- 
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ability faster as well as more photons ha ve been collected th an for the other sources. 
Current limits on Mqg > 3 x 10 17 GeV ( Albert et al.|2007i| are alre ady an order of 
magnitude better than limits obtained from GRBs (El lis et al.|2003 2006| l. It is im- 
portant to point out that if a dispersion of arrival times from flares is detected, source 
intrinsic spectral variations can be excluded once consistent time lags are detected 
from a sample of flares from different sources at different red-shifts. Ultimately, the 
inferred LIV characteristics should be consistent with a suppression of EBL absorp- 
tion processes. 



5 Concluding remarks on the future perspectives of 
the field 

The observational field of VHE gamma-ray astrophysics has been successfully driven 
in the last years by ground based imaging air Cherenkov telescopes. The "high en- 
ergy frontier" of Astrophysics will be expanded by the results expected from the 
Fermi satellite. The all-sky sensitivity of the Fermi Large- Area-Telescope (LAT) will 
most likely lead to many interesting new discoveries and surprises. The inter-play 
of the Fermi-LAT detections from space and follow-up observations from ground 
will provide mutual benefits for the communities as well as certainly answer many 
questions as well as stimulate new ones. 

While space-based observations of gamma-rays will be for a long time limited to the 
results obtained with Fermi (until maybe pair conversion telescopes will be deployed 
on the surface of the moon), ground based observatories will continue to be devel- 
oped and constructed during the next decade(s). To my knowledge, besides Fermi, 
no further space-based gamma-ray telescope, operating in the GeV energy range, is 
planned. Consequently, ground based gamma-ray detection techniques will be the 
only available experimental approach to observe high energy gamma-rays after the 
termination of the Fermi mission. 

The extensions of H.E.S.S. and MAGIC into phase II until 2009 will lower the en- 
ergy threshold and improve existing sensitivity moderately. The next generation of 
ground based installations w ill become fully operationa l at the end o f the next decade 



as envisaged in the AGIS 16 ( |Krawczynski et~aT|2007| and CTA 17 ( |Hermann||2007| ) 



projects. These installations aim at an improvement in sensitivity by a factor of 10 
and a widened reach in energy. At that time, th e new ground based instruments 



including larger (>km 2 ) installations like TenTen ( Rowell et al.|2008| will explore 
new energy windows above w 10 GeV as well as above 10 TeV 
However, there still remains a "blind spot": the transient sky above 10 GeV will nei- 
ther be explored very well with Fermi (limited photon rate) nor with ground based 
instruments of the current generation (limited field of view and energy threshold). 



New installations like HAWC (Sinnis & Hawc Collaboration 2005)1 will improve the 



situation in the energy range above a few 100 GeV. The lesson from the detection of 
fast transient events as seen from Cyg X-l and PKS 2155-304 however is that we are 



I6 Advanced gamma ray imaging system 
17 Cherenkov telescope array 
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very likely missing the fast variability of XRBs and AGN. Proper coverage of these 
objects requires instruments with a large field of view and an energy threshold well 
below 100 GeV. 
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